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We model and evaluate the Peltier and Seebeck effects in magnetic multilayer nanostructures 
by a finite-element theory of thermoelectric properties. We present analytical expressions for the 
thermopower and the current-induced temperature changes due to Peltier cooling/heating. The 
thermopower of a magnetic element is in general spin-polarized, leading to spin-heat coupling effects. 
Thermoelectric effects in spin valves depend on the relative alignment of the magnetization directions 
and are sensitive to spin-flip scattering as well as inelastic collisions in the normal metal spacer. 
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I. INTRODUCTION 



The Peltier effect refers to the conversion of an elec- 
tric voltage into a temperature difference (that can be 
used for refrigeration), while the Seebeck effect refers to 
the inverse process, the generation of an electric field 
by a temperature gradient. 1 Renewed interest in ther- 
moelectric properties is motivated in part by the im- 
proved performance of nanometer-scale structures 
Thin-film thermoelectric coolers can provide cheap and 
fast spot-cooling in micro- and nanoelectronic circuits 
and devices£>£ Strongly enhanced thermopower in quan- 
tum point contacts with widths approaching the Fermi 
wavelength can be used for sensitive and local electron 
thermometry.^* 9 - In ferromagnets and heterostructures in- 
volving magnetic elements, the effect of the magnetiza- 
tion (spin) degree of freedom on thermoelectric transport 
has to be taken into account . 10 ' 11 ! 12 The giant magneto- 
thermoelectric power in multilayered nanopillars 1 13 i 14 
thermally excited spin-currents in metals with embed- 
ded ferromagnetic cluster o 15 ' 16 and thermal spin-transfer 
torque in spin-valve devices^ 7 - are examples of spin- 
dependent thermoelectric phenomena on a nanometer 
scale. 

Recently a large Peltier effect was discovered in tran- 
sition metal multilayered nanopillars by Fukushima et 
a^ 18 i 19 The temperature and energy dissipation as a 
function of an applied current were monitored using the 
temperature-dependent electrical resistance. In asym- 
metric structures the parabolic dependence of the re- 
sistance arising from current-induced Joule heating was 
found to be modified by a superimposed linear (Peltier) 
term that shifts the minimum resistance to a finite value 
of the current that could be positive or negative, de- 
pending on the combination of materials. Gravier et 
aL— used a model of diffuse thermoelectric transport 
in (non-magnetic) metallic heterostructures to compute 
the Peltier effect. These calculations had to be carried 
out numerically and magnetism was not taken into ac- 
count. The sample cross-sections were used as fitting 
parameters that appeared to be too small compared to 
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FIG. 1: Schematic of a spin- valve structure connected to 
reservoirs at different chemical potentials and/or temper- 
atures ar e the Fermi-Dirac distribution functions). 
Charge and heat flows through the normal metal spacer (N) 
are functions of the angle between the magnetization direc- 
tions (9) in the ferromagnets. 



the actual sample sizes. 18 ' 19 This discrepancy was at- 
tributed to the neglect of interface scattering. Katayama- 
Yoshida et al£^ interpreted the perceived enhancement 
of the cooling power as a contribution from an adiabatic 
spin-entropy expansion term (fcs/e)ln2. Dubi and Di 
Ventral studied the Seebeck effect in single level quan- 
tum dots with ferromagnetic contacts. 

Enhancing the performance of solid state cooling 
elements remains a challenge both for theory and 
experiment.— i2iMi5iSiI Fukushima et al l 18 ' 19 suggested 
that the Peltier effect in transition metal nanostructures 
could be useful for cooling magnetoelectronic devices. 
In order to assess this idea, the material dependence 
of the Peltier effect in magnetic nanostructures has to 
be understood. In this paper we investigate the Peltier 
effect in magnetic heterostructures theoretically, taking 
into account spin-dependent interface and bulk scattering 
by means of an extended finite element (circuit) theory 
of transport- that is a generalization of magnetoelec- 
tronic circuit theory. 23 ' 24 ! 25 Such a theory is also suit- 
able to study the magneto-thermoelectric power in mag- 
netic multilayers in which transport is normal to the 
interfaces! 13 ' 14 

In the following we explain the method, initially dis- 
regarding bulk impurity scattering, which is justified in 
the thin layer limit when interface scattering is domi- 
nant. We start with a simple normal metal N\\N\N2 
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trilayer structure in which metal N is sandwiched be- 
tween two electron reservoirs consisting of different met- 
als Nin). We then extend the methodology by mak- 
ing first one and then both reservoirs magnetic. In the 
Fi\N\F2 (ferromagnet|normal-metal|ferromagnet) spin- 
valve sketched in Fig. [IJ the relative orientation of the 
magnetization directions as well as the strength of the in- 
elastic collisions become important parameters. Finally 
we show how bulk impurity scattering (normal as well as 
spin-flip) can be introduced into the formalism, discuss 
the relevance of our results for experiments and finish 
with a summary and conclusions. 

II. CURRENT-INDUCED ELECTRON 
COOLING AND HEATING 

A resistor network theory is an efficient way to describe 
the electric transport properties of magnetic heterostruc- 
tures. Used to model the giant magnetoresistance effect 
in terms of the two-channel series resistor model ) 26 ! 27 it 
has been generalized to include the spin transfer torque 
in non-collinear magnetization configurations.— A resis- 
tor model for thermoelectric effects was introduced by 
MacDonald 2 ^ to understand the effects of different types 
of impurities in homogeneous bulk metals. 

In the present section we develop a generalized (spin- 
less) circuit theory to describe the effects of temperature 
and voltage bias on normal metal heterostructures. Our 
starting point is a (non-unique) definition of the circuit 
or device topology by partitioning it into reservoirs, re- 
sistors and nodes. Discrete resistive elements are inter- 
faces, potential barriers or constrictions that limit the 
transport. We require that these resistors can be repre- 
sented by purely elastic scattering processes. The resis- 
tors are separated by nodes, in which electrons can be de- 
scribed by semiclassical distribution functions, in node 
i. When inelastic electron-electron or electron-phonon in- 
teractions in the nodes are strong enough, the fi become 
Fermi-Dirac distributions parametrized by temperatures 
Tj and chemical potentials /i^. The charge and heat cur- 
rents through a given resistor, denoted by I = J dei(e) 
and eQ = J de(e — /J,)t(e) = eI E — fil, respectively, where 
fj, is the global chemical potential in equilibrium, are de- 
termined via the spectral current density 

i(e) = G(e)[f L (e)-f R (e)}, (I) 

where G(e) is an energy-dependent (spectral) conduc- 
tance between two neighboring nodes with distribution 
functions Jl(R)- According to Landauer-Biittiker scat- 
tering theory the conductance 

G ( £ ) = x £ [Smn - r ™» {e)rmn (e)] ■ (2) 

mn 

depends on the energy-dependent reflection amplitudes 
r mn (e) at the resistor for electrons that are incident 
from a node or reservoir. 2e 2 /h ~ 1/13 kfl is the 



quantum of conductance. For structure elements such 
as interfaces which are not overly complicated, r mn (e) 
can be calculated from microscopic, first-principles 
calculations , 17 i 29 i 30 Thermoelectric effects arise from the 
energy dependence of G(e). The circuit theory approach 
requires that scattering in the resistive elements is elas- 
tic. G(e) may in principle be bias dependent, which 
becomes important for tunnel junctions. Here we con- 
centrate on abrupt intermetallic interfaces with bias- 
independent spectral conductances. The interface resis- 
tance of transparent interfaces in a diffuse environment is 
affected by a correction caused by the drift of the distri- 
bution function. 25 ! 31 The bare spectral resistance G(e)~ 1 
is then substituted by G(e) _1 - [GW £ (e) _1 + Gjv K (e) _1 ], 
where G Ni (e) = 2e 2 M Ni (e)/h are the Sharvin conduc- 
tances of the metals that form the interface and are 
the total numbers of single-spin transport modes. 

Let us now consider a non-equilibrium steady state 
in a simple 2Vi|JV|iV2 normal metal structure in which 
the chemical potential /Uip) = M — eVi(2) and temper- 
ature T\(2) of the nodes deviate from their equilibrium 
values (/j,,T). Node N is assumed to be fully ther- 
malized with a distribution described by (eV/v,T/v) that 
still have to be determined. To lowest order in the ap- 
plied thermoelectric fields, we use the expansions /$ ~ 
fo + If " V) + §?( T * " T ) in the following, where /„ 
is the equilibrium Fermi-Dirac distribution function. It is 
possible to proceed and compute thermoelectric proper- 
ties for arbitrary energy dependences of G(e). However, 
results become much simpler upon using the Sommerfeld 
expansion^ 2 - in ksT/eF- Provided the conductances do 
not vary too rapidly near the Fermi energy (to be precise 
d t G\ tF ^ oo and G(e)e- £ / fcsT | e>>eF -> 0), the following 
approximation is applied 

r dc cg{ ° ~ G{ ° )sn '° + Y {ksT)2 

x [„(„ - l)C- 2 G(0 + 2n£ n -%G(0 + ra|G(0] ?=0 

(3) 

where £ is the energy variable relative to the Fermi en- 
ergy. Applying the Sommerfeld expansion to the expres- 
sions for the charge and heat currents in terms of the 
spectral current density, Eq.fTJ), results in expressions 
for the charge and heat currents into the normal node 
through junctions 1 and 2 , valid for ^(k B T) 2 d 2 G\ eF <C 

G\e F - 

( qZ ) = Gl(2) ( ~ s m T Sl <S T ) ( Tn - S ) ' 

(4) 

where C = (k B /e) 2 ir 2 /3 ~ 2.45 x 1(T 8 V 2 K- 2 is the 
Lorenz number. Gi = Gi (cf) are the conductances and 
Si = —eCoTd t hiGi\ tF (Mott's formula) are the Seebeck 
coefficients or thermopowers at the zero-temperature 
chemical potential, which for metals is just the Fermi 
energy ep ■ In bulk materials the thermopower can be 
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positive or negative and even change sign as a function 
of the temperature^ 

The chemical potential shift eV/v and temperature Tjv 
of the central normal node and thus the thermal distri- 
bution function /jv are determined by conservation laws: 
the charge and energy flows of the electrons are con- 
served, I1+I2 — and Q1+Q2 = 0. The latter is affected 
in principle by the phonon heat conduction through the 
contacts (see Appendix A) but is disregarded here since 
thermal transport in good metals is dominated by the 
conduction electrons^ Using Eq. (01 and the conserva- 
tion laws we find to lowest order in the charge current / 
that the electron temperature Tjv of the normal island is 
modified from the zero-charge current value To as 



Tw — T 



(iii - n 2 )/ 

Kl + «2 



(5) 



The electron cooling or heating of asymmetric structures 
(Iii 7^ II2) by the applied current is the Peltier effect. 

n, 



1 i(2) — Si(2)T are the Peltier coefficients and 
K l(2) = A7 1 G 1 ( 2 )(1 



5*1(2) /A) 



(6) 



are the thermal conductances of the resistive elements. 
Large thermopowers violate the simple proportionality 
between electrical and heat conductance, K = CqTG, 
the Wiedemann-Franz law. The expression for the zero- 
current temperature in the central node 



To 



KlTl + K 2 T 2 
Kl + K2 



(7) 



follows from energy conservation. The charge current I = 
Iav+Iat is excited by a voltage difference AV = V2—V1 
as well as the temperature bias AT = T2 — Ti . The total 
electric charge and heat currents are then given by, 



I = G(AV + SAT) 
Q = -III - RAT, 
where the total conductance becomes 
GiG 2 (ki + k 2 ) 



G = 



G1K1 + (Gj + G 2 )K12 + G 2 k 2 



with 



_~ Gi Go / Si S? 
/sia = 2£ a T 1 [ 1 - ' 



Gi + G2 V A 



(8) 
(9) 

(10) 
(11) 



The total conductance/resistance (R — 1/G) violates the 
series resistor rule, R ^ R1 + R2 (with Ri — 1/Gj), which 
is only recovered when either Si — S2 or when Sf <C 
Cq; the latter is always the case at low temperatures. 
On the other hand, the total thermopower S (or Peltier 
coefficient Pi = ST) and thermal conductance R do obey 
simple sum rules 



(12) 
(13) 
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When in series, two thermal or electrical resistances are 
additive, as expressed in Eq. (fT3|) or R = Ri + i? 2 which 
is valid in the limit Sf <C Cq. In contrast, the ther- 
mopower, Eq. (j!2p . depends on the spatial distribution 
of the scattering objects rather than its integral (the ther- 
mopower in a homogeneous bulk metal does not depend 
on its length). Equation (fT2|) holds not only for the spa- 
tially distributed scatterers considered here, but can also 
describe the relative contributions of different types of 
scatterers to the thermopower in bulk materials.^ 

The second law of thermodynamics (a non-negative en- 
tropy production) requires k > 0, which in the Som- 
merfeld approximation leads to 3 -^ |5 ma;I ;| = -\/A — 
157/iV/K. Defining the thermoelectric figure of merit 
ZT = \SAT N /AV\, Eq. ([5]) for the temperature change 
ATjy = Tjv — T induced by an applied voltage results in 



ZT 



S(Si — S2) 



A - S 2 



(14) 



assuming Gi = G 2 and (Si — 5 2 ) Si(2) ^ A, so that 
l^maxl = a/A corresponds to the maximum efficiency 
ZT = 00. Even for the best thin film thermoelectric 
materials 3 - ZT ;$ 3, and for most metallic structures 
S 2 <C Cq. Quantum point contacts, however, have larger 
thermopowers due to size quantization^^ so S 2 can be 
comparable to Cq and the predicted effects should be ob- 
servable in nanoscale structures. 

Equation ([5]) holds for low current densities, for which 
SV/2 < n < 2C T 2 /8V, where 5V is the voltage drop 
over a single contact. Non-linear heating by applied 
currents can be included my means of the quadratic 
term in the expansions of the distribution functions, i.e. 

fi w /o+H (/i,-ri+f(T,-r)+0(ft-rf 2 +.., and 

using the approximation d£ (^^j £G(£) ~ G(0). 

When LT <g; 6V <C CqT 2 /H, expressions for the nonlinear 
currents reduce to 



Gi 



h(2) 

Ql(2) 

Vn - Vi(2) 
Tv - Ti( 2 ) 



(2) 



1 

— IIl(2) — ^1(2) — 



5*1(2) 
ATi( 2 ) 



(15) 



where V m = {V N + V m )/2 and T 1(2) = (T N + T 1(2) )/2. 
For reservoir temperatures Ti — T2 — T and in the ab- 
sence of the thermopowers 5*1(2) = 0, or when Joule heat- 
ing dominates, particle and energy current conservation 
requires (Gi + G%)Vn — G{V\ + G2V2 and 



T 



JY 



(AV) 2 GiG 2 
A (Gi + G 2 ) 2 



(16) 



so that we recover the result for the maximum ampli- 
tude of the electron temperature profile in the middle 
of a diffusive bulk wire (with the conductance Gi = 
G 2 = G) due to heating by inelastic electron-electron 
collisions! 36 ! 37 By taking into account the thermopowers 
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of the junctions, but in the limit S 2 <C Co, the change in 
the electron temperature in the island is found as 



AT N = 



(n x -n 2 )7 



RI 2 



C T(G 1 +G 2 ) 2£ T(Gi + G 2 ) 
/ AT\ 2 AT(AV) 2 S(AV) S 



\T J ' T £ T 2 ' £2 T 3 ' 



(17) 



in which the electric current Eq. ([5]) is driven by a voltage 
and/or temperature bias. Equations and (fTT|) yield 
the same result for a small temperature increase due to 
Joule heating when IT.! = n 2 and Ti = T 2 . 

Following Fukushima et a L 18 > 19 we use Eq. (fTT]) to 
derive an expression for the Peltier coefficient in terms of 
a critical current I p at which heating and cooling cancel 
each other. Assuming that the resistance scales linearly 
with the electron temperature in the node, AT\r = 
leads to 



RI 2 

(n 2 - n x )7 p = — * 



(18) 



The factor 1/2 on the right hand side implies that the cur- 
rent heating is only half as large as considered by Refs. 
[Hand[ll, whereas the "cooling power" RI p — 2(II 2 — IIi) 
(in units of mV) is twice as large. This discrepancy can 
be explained as follows. In our model, the energy is dis- 
sipated in the nodes and reservoirs of the device, not 
at the sharp interfaces. We monitor the temperature 
change in the normal metal node which is assumed to 
be effectively thermalized. Half of the generated heat 
is dissipated in the reservoirs that by definition do not 
contribute to the resistance change. The expressions of 
Fukushima et a/ . 18 ' 19 can be recovered by treating the 
highly resistive junctions in their samples as bulk ma- 
terial in which heat is generated and contributes to its 



temperature and resistance rise (see Section IVII 



III. PELTIER AND SEEBECK EFFECTS IN 
THE PRESENCE OF A SINGLE 
FERROMAGNETIC ELEMENT 

The thermoelectric transport Eq. (|4|) can be gener- 
alized to include the spin degree of freedom. For spin- 
dependent thermoelectric transport through an F\N in- 
terface the spin-polarized electric charge and heat cur- 
rents read 



jTU) 

Q 



TU) 



= G TU) 



1 STU) 

sni)T -c q t 



TU) 



JV 



TU) 

N 



U(i) 

ifu) 

(19) 

where the spin-dependence of the conductance G^^\ 
thermopower S^ = —eCoTd t In G^^\ eF , heat current, 
and temperature is expressed by the superscript j (J.) 
for majority (minority) spin electrons. V s = V~ — 
is the particle spin accumulation. Referring to the dis- 
cussion below we conjecture the existence of a heat spin 



accumulation T s = — T\ i.e., a temperature imbal- 
ance for majority and minority electrons, when thermal- 
ization is weak. We also define the total thermopower 
S = — (AV/AT)j = o of an interface between a normal 
metal and a ferromagnet as 



S = —eCoT 



d e G 
~G~ 



f G^S ^ +G l S l 
\ GT + Gi 



(20) 



This thermopower is observable when the interface is part 
of a (hetero) Sharvin point contact in direct contact with 
large reservoirs that prevent build-up of a spin accumula- 
tion. In a diffusive environment, however, the local spin 
accumulation should be taken into account, as described 
in the following. The spin-polarization of the interface 
thermopower is defined as 



Ps 



p'-p 
1 - P'P 



(21) 



where P and P' are the polarizations of the conductance 
G T ^) and its energy derivative d^G^^ respectively, both 
at the Fermi energy. Whereas \P\ < 1, \P'\ ^> 1 when 
d e G^ approaches —d e G^. Ps is also in principle un- 
bounded. Using 



(22) 



it follows that Ps ^ when the conductance polariza- 
tion is energy dependent. Spin-polarization of the ther- 
mopower of ferromagnetic material o 38 i 39 ' 40 has been in- 
voked to, e.g., explain the giant magneto-thermoelectric 
effect of magnetic multilayers^ For a few combi- 
nations of materials, the interface thermopower and 
its spin polarization are known from first principles 
calculations i 17 ' 30 

Consider now an F\\N\N2 pillar with one magnetic 
contact. Conservation of charge, spin and energy cur- 
rents implies the Kirchhoff rules I\ a + I 2ce = and 
J2 a Qi + Q2 — 0, where a =f (J.). The individual spin 
currents are separately conserved since we disregard spin- 
flip scattering in the normal metal spacer when the length 
of the metal does not exceed its spin diffusion length. 
In contrast, the heat spin accumulation, i.e. the tem- 
perature difference between the two spin species on the 
central island, is assumed to vanish by strong inelastic 
scattering, which is likely for temperatures which are not 
too low and/or metals which are not too clean. In this 
regime the electron temperature on the island becomes 



T; 



N 



We may call 



7i 



(71 5i - S 2 )I 
A>(Gi + G 2 )- 

G 2 / Gi 
(1-P 2 ) + G 2 /G! • 



To 



fl - P'P 



(23) 



(24) 



a " spin-entropy factor" , because it reflects the spin- 
polarization of the entropy flow per unit of the electric 
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current (the thermopower), S' 7^ S^ (P' 7^ P). In the 
limit S* a) < L (and therefore w LqTG^^) the 

temperature To (Eq. (JT])) is not affected by the mag- 
netism. When P' ^ P the Peltier cooling (heating) does 
not vanish even when Si — S2. The thermopower spin 
polarization, P$, can enhance or suppress the Peltier ef- 
fect depending on the spin polarization P and the relative 
amplitude of the conductances G2/G1. 71 can become 
large when \P'\ ^> 1 (an example is P' at a disordered 
Cr|Fe interface, see Table. I). 

The total thermopower can be expressed in terms of 
the properties of its constituent elements, in the limit 
S^Ci) <C \/2o an d for strongly thermalized electrons, as 



S 
G 



71 Si 
Gi 



G 2 ' 



(25) 



Therefore, when a spin accumulation is excited in the 
proximate normal metal, the magnetic junction con- 
tributes to the thermopower not by the Seebeck coeffi- 
cient of the point contact Si, but by the product with 
the spin-entropy factor 71 Si. 



IV. MAGNETO-PELTIER AND 
MAGNETOTHERMOPOWER IN SPIN VALVES 



We proceed to the study of thermoelectric effects in 
asymmetric Fi (mi) |iV|T2 ( m 2) spin valves (see Fig. [1]) 
for arbitrary relative orientations of the magnetizations, 
mi • rri2 = cos 6. The electron distributions in the 
nodes and reservoirs are now 2x2 matrices in spin space 
that can be expanded into scalar and vector components 

fF(N) = f F(N) 1 + & s F {N)f F(N)^ . . g ^ 

of Pauli matrices and 1 the 2x2 unit matrix. The unit 
vector of the spin quantization axis s F is parallel to the 
magnetization of the ferromagnet, whereas s N can point 
in any direction. In linear response, the 2x2 spectral 
current in spin space across a ferromagnet-normal metal 
junction at energy e in the absence of spin-flip and inelas- 
tic interface scattering is given as a spectral Landauer- 
Biittiker-like expressio n 23 i 24 i 25 



i N \F (e) = Yl G ° W ^ W " ^ 

a/3 



,0 



(26) 



where = (1 ± d"-m)/2 are projection matrices 

in which the unit vector m = s F denotes the mag- 
netization direction of the ferromagnet. The conduc- 
tance tensor elements read G Q/3 (e) = (e 2 /h) z^ nm [8 m n — 
r nm ( e ) ( r nm ( e ))*] m terms of the energy-dependent re- 
flection coefficients r" m (e) for majority and minority 
spins at the N\F interface. Its diagonal elements are the 
conventional spin-dependent conductances that govern, 
e.g., the giant magnetoresistance, whereas the complex 
non-diagonal elements, the so-called spin-mxing conduc- 
tances, parameterize the transverse spin currents that 
are absorbed by the ferromagnet and give rise to torques 



on the magnetization. The total charge-spin and heat 
matrix currents are defined as I = J de ?(e) and eQ = 
J de(e — jti)i(e) = eP — fit, respectively, where /i is the 
equilibrium chemical potential and P the energy cur- 
rent. In the following we assume that both spin com- 
ponents /Tw = f c ± f s of the diagonalized matrix dis- 
tribution functions f F ( N > may be described by thermal- 
equilibrium Fermi-Dirac distribution functions with spin- 
dependent chemical potentials and temperatures. The 
Sommerfeld expansion can then be employed to derive 
expressions for the transport currents as a function of 
applied voltage or temperature gradients in terms of the 
conductance tensor G af3 and its energy derivative G™' 3 at 
the Fermi energy^ The total charge, spin and heat cur- 
rents read I c = Tr[I], I s = Tr[oT] and Q c = Tr[Q], re- 
spectively, where the trace is over spin indices. 

The charge and energy conservation laws read I c i + 
I C 2 = and Qd + Q C 2 = 0. Moreover, in the ab- 
sence of spin-flip scattering in the normal node, the to- 
tal spin angular momentum current is conserved as well, 
i.e., I s i + I S 2 = 0. These Kirchhoff Laws close the sys- 
tem of transport equations in the strongly thermalized 
regime. In what we call the weakly thermalized regime, 
the distributions for each spin species are thermalized 
separately, but the energy exchange between the spin 
subsystems is disregarded, which is a realistic scenario 
at low temperatures^ 2 " In this limit a spin temperature 
vector on the central island exists, T s ^ 0, and we re- 
quire Q s i + Q S 2 ~ where Q s = Tr[<xQ], which means 
that energy is conserved for each spin channel separately. 
It is worth while to compare the thermoelectric trans- 
port properties such as the total conductance and ther- 
mopower of the spin-valve structure in the different in- 
teracting regimes. 

In the strongly thermalized regime for a symmetric 
spin- valve, G°f — G^ and G"f — G"^ , the tempera- 
ture of the normal metal island Tv = To is not affected 
by electric current. The total electric current readsii 

h = f (AV + SAT)-^ 'T^ilfo ( PAV + P ' SAT ) ' 
2 2 r)R + tan 0/2 

(27) 

Here t)r = 2ReG^/G > 0, where G^ is the complex 
spin mixing conductance^ For most metallic contacts 
r]i = 2 Im G^-l /G, is small (rjj <C Vr)^ an d is disregarded 
in the analytical results. However, since 77/ is not small 
for Cr|Fe and Cr|Co interfaces, see Table I, it is included 
in the numerical results for these junctions. The angular 
magneto-resistance for AT = as measured by Urazhdin 
et al£& is well described by circuit theory4£ The thermo- 
electric transport properties of the spin-valve structure 
differ significantly in the different interacting regimes. 
In the Sommerfeld approximation, the spin-mixing ther- 
mopower STi = _ e£ TGP/G^ and the dimensionlcss 
mixing parameter = 2Re(Im)GP/G e enter expres- 
sions for the electric currents only when T^ ^ 0, i.e. in 
the weakly thermalized regimeJ^ 
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0.48 


Cu|Co(lll)* 


3.42 


-15 


68 


77 


17 0.64 


-0.073 


0.13 


-0.05 


Cr|Au(001) 
or|/vu^uui. ) 


0.36 
u.o / 


7 

n 
U 










- 
- 


— 


— 


— 


Cr|Fe(001) 
Cr Fe(OOl)* 


0.88 
0.94 


22 
7 


-74 
-53 


-40 
-190 


48 4.23 
-9500 3.25 


1.38 
0.43 


-4.27 
-0.48 


-1.38 
9.39 


Cr|Co(001) 
Cr Co(OOl)* 


0.56 


62 


-62 


-111 


-160 3.03 


-0.59 


-2.86 


-3.46 


0.71 


23 


-23 


-95 


-92 2.92 


-1.79 


-0.86 


3.21 



TABLE I: Thermoelectric interface parameters calculated at the Fermi energy for a number of almost lattice-matched interfaces 
including Schep's drift correction.— With the exception of the Cr|Co interface for which Co is assumed to be bcc, Cr and Fe are 
bcc while Cu, Co and Au are assumed to be fee. The asterisk * indicates a dirty interface modeled in a 10x10 lateral supercell 
with two layers of 50%-50% alloy. 



In the presence of a temperature bias AT with an open 
electric circuit (J c = 0), the induced thermoelectric volt- 
age AV^ is described by the angular magnetothermopower 
(MTP) S(9) = - (AV/AT) Ic=0 which in the strongly 
thermalized regime readsii 



S(0) _ m + (1 - PP') tan 2 9/2 
?te + (l-P 2 )tan 2 6>/2 ' 



S 



MTP: 



S(n) - 5(0) 
5 



-PPs 
1 + PPs' 



(28) 



(29) 



The MTP is finite when the interface thermopower is spin 
polarized, P$ ^ 0. When PP' > 1 (which also requires 
P s < 0) one finds an angle O = 2 tan" 1 y/r) R /(PP' - 1) 
where the thermoelectric power can change sign, i.e. 
a spin-valve with non-collinear magnetic configuration 
can display a transition from electron-like to hole-like 
transport. In a spin- valve the thermally induced charge- 
currents that enter the normal node from both sides can 
be made to cancel such that the net thermoelectric volt- 
age vanishes, 5 = 0. Note that not only the individual 
spin currents but also the charge current depend on the 
effective spin polarizations. The MTP vanishes in the 
half-metallic limit P = P 1 = ±1. 

In the weakly thermalized regime on the other hand, 
the magnetothermopower Eq. (|2"9")l is twice as large, 
whereas the magnetoresistance and 5(0) do not change 
(in the limit 5 2 <C Co). S(n) is enhanced in this case 
because the heat spin accumulation facilitates the ther- 
moelectric voltage build-up. In this regime an angular 
magnetothermopower is found even when P$ = pro- 
vided that rj' R = 2Red e G^ /<9 e G ^ rm, , which is de- 
stroyed by full thermalization. 

The Onsager-Kelvin relation between the total Seebeck 
and Peltier coefficients in spin valves, i.e. 11(0) = S(0)T 

in which ft = (— Q c /Ic) , is found to hold in both 

V _ / at=o 
thermalization regimes. 



For a quantitative analysis we need to know more 
about the thermoelectric parameters for interfaces. In 
the absence of experimental estimates, we calculated the 
parameters from first-principles within the framework of 
density functional theory for a number of interfaces which 
figure prominently in the field of magnetoelectronicsJ^ 
The values are given in Table fl] For an A|B interface, 
the calculation proceeds as follows. 29 ' 46 Self-consistent 
density functional theory calculations are first performed 
separately for bulk A and B materials. These calculations 
yield bulk charge- and spin- densities and potentials, and 
the corresponding Fermi energies. A self-consistent inter- 
face calculation is next performed subject to the poten- 
tials (and densities) far from the interface being equal to 
their bulk values, up to a constant which is adjusted so as 
to equalize the Fermi energies.— The interface breaks the 
lattice periodicity perpendicular to the interface leaving 
only two-dimensional periodicity parallel to the interface 
that is characterized by the two-dimensional Bloch vec- 
tor k|| . The electronic structure of the localized pertur- 
bation formed by the interface is handled using a Green's 
function method, a so-called "Surface Green's Function" . 
The rank of the matrix of the perturbation is made finite 
and minimized by making use of the translational sym- 
metry parallel to the interface and using a maximally- 
localized basis of tight- binding (TB) muffin-tin orbitals 
(MTOs). 48 ' 49 To calculate the scattering matrix S 



S(k 



_fr(k h e) i'(k||,e) 
t(k|,,e) r'(k| h e) 



(30) 



at real energies (at or close to the Fermi energy in the 
context of transport), we use a wave- function-matching 
scheme due to Andc 50 which involves the calculation of 
individual scattering states far from the interface. The 
rank of the reflection and transmission matrices r, r' , t, t 
is determined by the number of Bloch states at a given 
energy e and transverse wave- vector k|| . The mini- 
mal TB-MTO basis is very efficient and makes it possi- 
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FIG. 2: Magnetothermoelectric power as a function of relative 
angle between the magnetizations in a Fe|Cr|Fe (001) spin 
valve with dirty interfaces. The MTP is significantly different 
for the strong (full line) and weak (dotted line) thermalization 
in the normal metal spacer. 



ble to model incommensurate lattices and various types 
of disorder using large lateral supercells i 29 ' 46 ' 51 Substi- 
tutional disorder where one or more layers of atoms 
form an alloy is conveniently treated by calculating 
the potentials self-consistently using a layer version 4 ^ 
of the coherent potential approximation^ and then dis- 
tributing at random the site potentials in lateral su- 
percells subject to maintenance of the appropriate layer 
concentrations^ The mixing conductance is most eas- 
ily calculated in terms of the reflection matricesj 43 i 53 We 
consider here interfaces in diffuse metallic systems which 
implies that we have to use a generalization of the Schep 
correctio n 25 ' 31 by replacing the bare G Q ^(e) _1 with 
G^ie)- 1 - [G N (e)- 1 +-5 Q/3 Gp(e)- 1 ], where G N (e) = 
e 2 M N (e)/h and Gg(e) = e 2 M§t(e)/h are the single-spin 
Sharvin conductances of the normal and ferromagnetic 
metals forming the interface. The thermopower and 
other generalized thermoelectric parameters are deter- 
mined by numerically differentiating the scattering ma- 
trix calculated as a function of the energy. Details of the 
numerical procedures will be given in a separate paper ^22, 

We use the data in Table I to compute the angular 
dependence of the thermoelectric properties of a few spin 
valves for illustrative purposes. Spin-flip scattering is 
disregarded here but will be discussed in later sections. 

We plot the angular MTP for a Fe|Cr|Fe (001) spin 
valve with dirty interfaces in Fig. [5] and for a Co|Cu|Co 
(110) spin valve with clean interfaces in Fig. [3] The an- 
gular MTP is enhanced in the weakly thermalized regime 
(shown by dotted lines) by up to a factor of two for the 
antiparallel configuration. Depending on P and Ps, the 
MTP can be of any sign, see Eq. PTJ|) . 

We now turn to the Peltier effect of asymmetric spin 
valves. In the strongly thermalized regime and for Ps t — 
Ps 2 = 0, the Peltier cooling retains the simple form 
for normal metal structures, Eq. (0, whereas the to- 




FIG. 3: Same as Fig. [2] but for a Co|Cu|Co (110) spin valve 
with clean interfaces. 



tal charge current is a complicated function of the mag- 
netic configuration of the system. A magneto-Peltier 
effect (MPE), i.e., a dependence of the cooling power 
on the magnetic configuration, is found when the ther- 
mopower is spin dependent. For a voltage-biased spin 
valve with thermal asymmetry Si =/= S 2 and P[ 7^ P 2 , 
but G x = G 2 = G, Pi = P 2 = P, and Sj {l) < y/Z^, 
we find for the temperature change of the normal metal 
spacer 



ATjv = 



2£ G 



where the spin-entropy factors 

-P 1(2) P) tan 2 9/2 + m 



7l(2) V>) 



(1- 



-P 2 )tan 2 0/2- 



- VR 



(31) 



(32) 



depend now on the magnetic configuration {P[/2) =h P)- 

7^ = 1 and (l - P[ (2) P^j I (1 - P 2 ) , respectively, for 
parallel and antiparallel configurations. The MPE should 
therefore be observable in R vs. I curves of spin valves 
during current-induced magnetization reversal. Accord- 
ing to Eq. (18]) with R(I = I p ) = R(I = 0) = Ro, the 
temperature change AT/v = Tn(6) — To corresponds to 
the cooling-power RqI p ~ 2£ (Gi + G2)TATn/I p such 
that 



R [1,(70-1,(0)] - 
G 1 P 2 {P{-P 1 )U 1 - 



G 2 p 1 (p^-p 2 )n 2 



G 1 (l-P 2 ) + G 2 (l-P 2 ) 



(33) 



This signal contains unique information on the spin- 
polarization of the thermopower. When thermalization 
is weak a MPE arises even when Pg 1 — Ps 2 = 
{P[(2) = ^1(2))- A sign change in the cooling-power is 
also expected. In the strongly thermalized regime this 
arises from different angular dependences of the spin- 
entropy factors. When the effective thermopowers are 
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FIG. 4: Magneto-Peltier cooling as a function of the angle 
between the magnetizations in the hypothetical bcc Co|Cr|Fe 
(001) spin valve. RqI p (8) is the cooling power which com- 
pensates local Joule heating. The signal Ro[I p (n) — I P (0)] is 
significantly different for strong (full line) and weak (dotted 
line) thermalization in the normal metal spacer. 



equal, 7 1 MP ((9 )S'i = ^ P {6 )S 2 , no Peltier cooling is ex- 
pected, RoI p {0q) — 0. 

In Figs.[4]and[5j we illustrate the magneto-Peltier cool- 
ing by computing the angular dependent cooling-power 
RqI p (9), at room temperature, for a hypothetical bcc 
Co|Cr|Fe (001) spin- valve structure with clean interfaces 
and for an asymmetric Co|Cu|Co (001) with one ideal and 
one disordered interface, using parameters from Table I. 
For comparison, we include the results in the weakly ther- 
malized regime, depicted in the figures by dotted lines. 
The dependence of the Peltier cooling on the magnetic 
configuration observed in Fig. Q]for the Co|Cr|Fe spin- 
valve is caused by the relatively large values of the inter- 
face thermopower parameters P and P' . The relatively 
weak magneto-Peltier signal, Ro[I p (n) — I p (0)] vanishes 
when Gi = G 2 and P\ — P2, for weakly thermalized 
electrons (dotted line) is strongly enhanced by the inter- 
spin energy exchange in the opposite limit. The MPE 
for the asymmetric Co|Cu|Co (001) structure, Fig. [5j 
displays similar effect but with smaller amplitudes. For 
the strongly thermalized electron case, the MPE caused 
by interface scattering, Rq[I p (tt) — I P (0)]> is of the or- 
der of 1 — lOmV, which is smaller than the experimental 
value o 18 ' 19 RqI p ~ 20 — 40mV. For a typical current 
density of 10 8 A/cm , we find a maximum temperature 
drop Aljv(fl) ~ 0.15K which is also too small to explain 
experiments. Note that in the above we have assumed 
that interface scattering is dominant. As explained in 
the next sections, we show that these numbers are in- 
creased by including scattering in the bulk. 

The magneto-Peltier effect vanishes for symmetric 
spin-valve structures. Introducing the asymmetry G\ 7^ 
G 2 , we obtain an angle-dependent temperature modula- 



R I p [mV] 




FIG. 5: As Fig. g]but for an asymmetric Co|Cu|Co (001) spin 
valve with one clean and one disordered interface. 

tion 

£ o AT N {0) = 

2{G 1 -G 2 )P (P' - P) SI C sin 2 § 
(d + G 2 ) 2 (l - P 2 ) VR + GiG 2 (l -P 2 - VR? sin 2 0' 

(34) 

even for equal thcrmopowers and spin polarizations at the 
interfaces. Analytical expressions for the Peltier cooling 
in the weak thermalization regime, in which the spin- 
mixing thermopower (f^n in Table I) becomes a rele- 
vant parameter, are much more complex. The computa- 
tion is straightforward, however, and is easily carried out 
when the necessity arises. 

V. SPIN-CONSERVING BULK IMPURITY 
SCATTERING 

In this section we discuss the contribution of bulk scat- 
tering for the case of wires with constant cross section 
Ao, first for non-magnetic metals and then for magnetic 
structures, both in the strongly thermalized limit. 

A normal metal pillar is a heterostructure 
K A \N A (L A )\N(L)\N B (L B )\K B , where N^U) de- 
notes a layer of material i (=A, B) with thickness 
that can be larger than the elastic mean-free-path due to 
disorder scattering. The length L of the central island N 
is so short that its (bulk) resistance can be disregarded. 
The external reservoirs TZa(b) & r e in thermal equilib- 
rium but at different temperatures and/or voltages. 
The spreading resistance at an abrupt opening can be 
accounted for by an effective length parameter^ The 
electron distribution functions in the disordered metal 
wires follow from the diffusion equation in the bulk and 
are connected at the interfaces by (quantum mechanical) 
boundary conditions.— The conserved particle/heat 
currents can be obtained from Eq. ^ by replacing 
-(^1(2) - V N ) by L A{B) V^ A(B) /e , T 1(2) - T N by 
La(b)VT A (b) and the interface conductance Eq. @ by 



the electric conductivity (Ja(B) — e2 -N~A(B)D a(b), where where the spin-entropy factors 7^(#) are 



Af and D are energy-dependent densities of states and 
diffusion constants of the bulk materials, respectively. 
Mott's formula, Sa(b) = — e £oTd e hiaA(B)\e F , holds 
for the diffusion thermopower which usually dominates 
at high temperatures ! 54 ' 55 In linear response charge 
and energy current conservations imply V 2 m = and 
V 2 T = 0. The chemical potential and temperature 
depend linearly on position except for jumps at the 
contacts which are governed by the interface parameters. 
The local chemical potential and temperature are then 
found by the charge and energy current conservations at 
the boundaries. 

As a function of the applied electric current we obtain 
the following expression for the temperature change AT/v 
on the normal metal island 



AT N = 
( Sa — Sb 



S\ — 02 . Sa — S2 , S\ — Sb \ GtotI 



\ GaGb G\Gi 



GaG2 



G\Gb 



Co ' 
(35) 



where Ga(b) = &A(B)Ao/ 'La(b) and Sa(b) are the bulk 
(Drude) conductances and thermopowers in the leads, 



G: 



2.A.B 



G: is the total scries conductance. 



The interface contribution to the Peltier cooling disap- 
pears when Si — S2 and G1/G2 — Ga/Gb- When Sa — 
Sb and a a = &b > Peltier cooling is possible for different 
lengths of the normal leads (La ^ Lb)^ Eq. ([33)1 can 
be simplified by introducing the lumped conductances 
G L = GaG\/(Ga + Gi) and G R = G 2 G B /{G 2 + G B ) 
as well as the thermopowers Sl and Sr with Sl/Gl = 
Sa/Ga + Sx/d and S R /G R = S 2 /G 2 + S B /G B for the 
left and right parts of the normal island. In terms of the 
new parameters we find 



AT, 



JV 



{S L ~ S R )I 
£o(Gl + Gr) 



(36) 



which, as expected, has the same form as Eq. ([5]) in the 
limit S 2 < C . 

Replacing the normal lead Na by a magnetic lead, 
say Fa, we find that the thermoelectric cooling obeys 
Eq. (|2"3")l after replacing the interface conductances and 
thermopowers Gi(2) and «Si(2) by Gl(r) and Slir), pro- 
vided that the spin polarizations in bulk layers and con- 
tacts are the same. A more complicated structure like 
Co|Au|Ti|Au can be shown to be equivalent to an F|N|N2 
pillar by a similar lumping of parameters. 

We now turn to the MPE, i.e. the dependence of the 
Peltier cooling on the magnetic configuration of a spin- 
valve structure, in the presence of bulk scattering. A 
simple analytical expression for the cooling-power (or the 
local Joule heating compensation current I p ) can be ob- 
tained when the spin polarizations of the bulk and inter- 
faces are equal: 



7f P 



G L [1 - P A P B - {Pa - Pb) P' a \ + G R (1 - Pf) 
G l {1-P\)+Gr{1-PI) 



An expression for ~{ R P is obtained by interchanging the 
indices L <-> R and A <-> B. Eqs. (|37l38p reduce to Eq. 
(|33f when bulk scattering is disregarded. 

A phonon (or magnon) thermal current can transfer 
momentum to the electrons in the presence of inelastic 
scattering which in turn generates an additional electric 
field and modifies the thermopower. For normal (as well 
as ferromagnetic) metals at sufficiently low temperatures 
a contribution of the phonon (and magnon) -drag effect 
may become significant ! 58 ' 59 The magnon-drag effect is 
likely to be suppressed strongly in heterostructures since 
magnons cannot escape the ferromagnets. Strong phonon 
scattering at interfaces will likewise reduce the phonon- 
drag effect in multilayers. A microscopic treatment of 
the phonon-drag effect in heterostructures is beyond the 
scope of the present paper, however. At elevated temper- 
atures, where the drag effect can be disregarded, Mott's 
formula holds approximately even in the presence of in- 
clastic scattering! 56 ' 57 



VI. SPIN-FLIP BULK IMPURITY 
SCATTERING 

Here we study the influence of spin-flip relax- 
ation on the Peltier and Seebeck effects in magnetic 
K A \F A {L A ) \N (L) \F b (L b ) \K b nano-pillars, where F A 
and Fb denote disordered ferromagnetic layers , with 
collinear magnetization directions. We assume that bulk 
impurity scattering is dominant so that interfaces may be 
disregarded. The charge and spin distribution functions 
in the ferromagnet, / c ( s ) = (/ T ± f^) /2, respectively, are 
then solutions of the spin diffusion equations that are 
continuous at the interfaces^ In the strongly thermal- 
ized regime, defining /i s = /i| — /ij and /i c = (/ij + /i|)/2 
as, respectively, the spin and charge chemical potentials, 
we find (see Appendix B for details) the following ther- 
moelectric spin diffusion equations in a ferromagnet 



V 2 /i s 



12 ' 



VVc = -P^t, 
Zl sf 

V 2 T 



(P'-P)S ^ 



C 







2/ 2 ' 



(39) 
(40) 

(41) 



R I P = 2( 7 f ^n L 



(37) 



Here l s f stands for the spin-flip diffusion length and P 
and P' are the spin polarizations of the bulk conductiv- 
ity and its energy derivative in the ferromagnet. These 
equations have to be solved with continuity boundary 
conditions at the interfaces. The expressions for the cur- 
rents are similar to Eq. (|19p after replacing temperature 
and voltage differences by gradients and conductances 
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by conductivities. Eq. (|41j) has to our knowledge not 
been given elsewhere, but is required by the conserva- 
tion of charge and energy currents. According to this 
equation the decay of the spin accumulation in the ferro- 
magnet provides a source or sink of heat currents (when 
P' 7^ P). We can understand this effect by the charge 
accumulation that is locally generated by spin flips in fer- 
romagnets, Eq. (|40p . Similarly, spin- flip scattering in the 
presence of spin polarization of thermopower modifies the 
distribution functions in a way that can be interpreted 
as a source or sink of heat as expressed in Eq. (|4T~j) . The 
spatial variation of /x c ( s )(x) and T(x), in a voltage biased 
spin- valve is sketched in Figs. [6] and [7] The local charge 
and spin chemical potentials in a spin-valve biased with 
a voltage difference (Fig. [6J do not depend much on the 
thermopower or the Peltier cooling. More interesting are 
the results in Fig. [7l illustrating the strong dependence 
of the local temperature on magnetization configuration, 
the strength of the spin-flip scattering and the spin po- 
larization of the thermopower. 

In the following we assume identical spin polarization 
and spin-flip diffusion length l s f in the magnetic leads 
and l^j S> L. Our results become simple when the only 
asymmetries of the pillar are a a ^ <*b , Sa ^ Sb- For 
parallel alignments of the magnetizations the Peltier cool- 
ing is equal to that of a normal metal structure with 
R o I P {0) « 2£ (G A + G B )TAT N /I p = 2(U A -U B ). How- 
ever, the magneto-Peltier signal in the presence of spin 
decay becomes 



Ro[I P {7r)-I p {0)] = - 



tanhA 4PP S G A U A -G B U B 
A 1 + PP S G A + Gi 



where A = La/I a 



(42) 

s j — L B /l^j is a measure of the spin- 
flip scattering in the ferromagnets. The magneto-Peltier 
signal decays with increasing A, e.g., on using a thicker 
magnetic leads. The magneto-Peltier signal vanishes 
when A » 1, and reduces to an expression equivalent 
to Eq. (|33[) in the opposite limit. Spin- flips in the nor- 
mal metal spacer (with thickness comparable to or longer 
than ifj) also reduce the magneto-Peltier signal. 

For the spin-valve structure in the presence of bulk 
spin-diffusion in the ferromagnets, we find the following 
results for the magnetoresistance MR as well as the mag- 
netothermopower MTP. For the parallel configuration, as 
expected, no spin-flip contribution to the total resistance 
and thermopower is obtained i.e., Rp — Ra + R B and 
Sp = (R A S A + R B S B )/(R A + R B ). However, for the 
anti-parallel configuration we find 



Rap — Rp = 



tanhA 4P 2 R A R B 



Sap — Sp — — 



A 1 
tanh A 



-P 2 Ra 
2PP S 



-R B 
RaRb 



A l + PPs (Ra + Rb)' 



:(S A 



(43) 

S B ) 
(44) 



The MTP is therefore proportional to the giant 
magneto-resistance, independent of the spin-flip scatter- 
ing strength A. In contrast to the MR (at A 3> 1) it is 
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FIG. 6: Illustration of the local charge and spin chemical 
potentials in a Fa|N|Fb spin valve biased by a voltage dif- 
ference, for both parallel (TT) an d anti-parallel (jj) align- 
ments of the magnetizations. Fb(P = 0.8, As = 0.1, Sb 

1 pV /K) has been chosen to have weak spin- flip scattering 
and thermopower compared with Fa (J 3 = 0.8, \a ~ 10, Sa = 
— 20/iV /K). The thin normal metal spacer is chosen to be 
highly conductive (pjv <S Pa(b) = 10/^ficm) and to have a 
thermopower equal to that of Fb- Note that spin accumula- 
tion is assumed to vanish at the two ends (at reservoirs). 
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FIG. 7: Illustration of local temperature distribution in the 
voltage biased Fa|N|Fb spin valve structure as in Fig. |S] 
for both parallel (jj) and anti-parallel (jj) alignments of 
the magnetizations. The different temperature profiles corre- 
spond to (a) P' < P (P s = -0.8), (b) P' > P (P s = -1.8). 
For the parallel alignment of the magnetizations the Peltier 
cooling is insensitive to spin-flips and reaches its maximum in 
the central node. 



not possible to model spin-flip scattering for magneto- 
thermoelectric effects by replacing Ra(b) by the resis- 
tance of the magnetically active region A -1 i?vum- 
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VII. RELEVANCE FOR EXPERIMENTS 

In the magnetic nanopillars considered by Fukushima 
et a h ls ^ 19 the magnetic or normal leads that connect the 
central spacer to the wide external reservoirs are so long 
that the bulk scattering is important: A 100 nm long 
cobalt wire has a resistance pcol ~ 6 f £1 m 2 at room 
temperature, which is larger than the interface resis- 
tance Rco\Cu ~ 0.25 iVLm 2 (see Table. I). The effec- 
tive thermopower Sl ~ Sco — —31 fj,Y /K from the bulk 
scattering is also larger than the interface thermopower 
•ScolCu ~ — 6 fiV /K. The interface contribution to the to- 
tal thermopower would become more important for high- 
resistance interfaces, such as tunneling barriers or point 
contacts, or structures with thinner layers. 

The lattice exchanges energy with the conduction elec- 
trons by inelastic electron-phonon interactions. In prin- 
ciple, there is a net heat current flowing between the elec- 
tron system and the lattice/substrate. In a steady state 
situation it is reasonable to assume that electron and 
lattice temperatures are identical and resistance changes 
reflect the electron temperature. When Peltier cool- 
ing and Joule electron heating compensate each other 
the temperature change vanishes. Estimating the non- 
linear electron heating in the island by CqTATn = 
Gt ]l 2 /2(G L + G R ) and using Eq. (23) in terms of the 
lumped conductances and thcrmopowers, the compensa- 
tion current I p or cooling-power RqI p for F.a|N|Nb struc- 
tures such as Co|Cu|Au nanopillars can be expressed as 

R I P «2( 7L n L - IL R ). (45) 

Such an expression holds as well for 
Co|Cu|Co|Au structures when the magneto- Peltier effect 
for the symmetric part Co|Cu|Co can be disregarded. 1 ^ 19 
The factor of 2 difference with Refs. Il8lll9l has been 
noted already above. We find below that including this 
factor leads to a better agreement of a simple model 
of bulk thermopowers with experiments. Also the too 
large pillar cross sections with which Gravier et aZ.— 
fitted their numerical results to the experiments can 
be traced back to this factor 2 in the Joule heating. 
Our model might not be appropriate for Fukushima's 
samples that contain a highly resistant, presumably 
oxide, layer over which much of the voltage drop occurs. 
Such a layer, when sufficiently thick, might be better 
described as a bulk resistor in which Joule heat is 
preferentially generated. At the compensation current 
I p , finite temperature variation profiles may persist since 
the Joule and Peltier sources are spatially separated. 
A simulation beyond our simple model might then be 
required for a quantitative description. 

Experimental values of RqI p for Co|Au nanopillars can 
be read off the figures published by different groups, 
amounting to (in mV) 19 (Ref. M), 23.0 (Ref. M), 
22.5 (Ref. HUli). These numbers agree well with the 
following results. For a finite length of the bulk lay- 
ers, i.e. Leo = 60 nm and Lau — 120 nm and 



taking into account the interface scattering, when disre- 
garding the spin polarization of the thermopower Pg we 
find RqI p — 19.5 mV . Here we also assumed Gr «Gi 
caused by an oxide layer on the non-magnetic side of the 
structure, Scu|Au — 0, ffc |Qu — —6 fiV /K and used the 
bulk parameters from Ref. [2(j|. RqI p ~ 2(IIco — lL4„) = 
19.6 mV indicates that the Peltier cooling is not sig- 
nificantly affected by interface scattering. A finite Ps 
can enlarge or reduce the above estimates. The spin- 
entropy coupling factor j L w (1 — P'P)/(1 — P 2 ) when 
the bulk and interface spin polarizations are the same. 
For Co wc took P = 0.44 (Ref. Hi). Conflicting val- 
ues P s = -0.18 < (S T = 0.75;) (Ref. and 
P s = +0.42 > (S T = -30fiV/K,Si = -12/xV/K) 
(Ref. HI) are found in the literature. According to Table 
I, Ps of the Co|Cu interface can also have either sign. The 
two values for Ps modify the above estimate to RqI p = 21 
and 17 mV, respectively, possibly favoring a Ps < when 
comparing with the observed values. 

The adiabatic spin-entropy expansion term 
(kBTln.2)I /e considered in Ref. [2l| is in our opin- 
ion an extrapolation of a concept from equilibrium 
thermodynamics that does not play a role in the current 
induced (non-equilibrium) Peltier cooling. 

We proceed by estimating the magnitude of 
the temperature drop that can be realized by 
the Peltier effect in the magnetic heterostructure 
Co(60nm) |Au(20nm) |Cr(120nm)^ At room tem- 
perature the bulk thermopowers of both Co and 
Cr are relatively large and have opposite signs 
(Scr = +21.56 /xV/K). The temperature drop 
AT Au « I P {S L - S R ) I {Co [G L + G R )) in the cen- 
tral island amounts to 4.8 K at I p = 10mA for a 
cross-section of 70 nm x 200 nm, at a current density 
of ~ 10 8 A / cm 2 , which is close to the maximum tem- 
perature drop in the temperature profiles computed in 
Ref. (we find a cooling-power RoI p — 30 mV, which 
is smaller than the observed value 41 mV, however). 
The temperature reduction per unit of electric current 
is sensitive to the thickness of the leads. For the thick 
magnetic layers ATa u ~ I P (Sco - Sbr) / £ (Gco + Go). 
Spin polarization of the thermopower in Co can modify 
the amount of the temperature reduction, up to 8% for 
\Ps\ = 0A. 

In Fukushima's experiments the leads connected 
to the external reservoirs are long compared to 
the spin-flip diffusion length. In that regime a 
magneto-Peltier effect should be small. Let us 
therefore consider a spin-valve structure such as 
Coooi (10 nm) |Au (20 nm) |Coooi (5 nm) , in which spin- 
flip scattering is less important. Due to the different 
lengths of the bulk Co layers the Peltier cooling does 
not vanish at this structure even for a parallel mag- 
netic configuration; recall for example Eq. (|34[) when 
G\ 7^ G2. For the parallel alignment of the magneti- 
zations with Pa = Pb one finds I^r) = ^' equivalent 
to a normal metal structure, whereas the spin-entropy 
coupling parameter differs for the anti-parallel magnetic 
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configuration when Pa — —Pb- Let us now consider 
a small G\ = 0.01G 2 e.g. caused by an oxide layer at 
the junction between the thick Co layer and the normal 
metal spacer, using data in Table I for the interface scat- 
tering (at room temperature), and adopting bulk val- 
ues Ps = —1.18 and Ps — 0.42, we find respectively 
the magneto-Peltier signals Ro[I p (n) — I p (0)] = —1.59 
and +6.7 mV, which should be experimentally observ- 
able. Replacing the bulk parameters of the thicker Co 
layer by p Fe = 9.7 x l(T 8 ftni and S Fe = +20^V/K, 
the Peltier cooling is increased and the magneto-Peltier 
signals read R [I p (tt) - I p (0)] = +3.2 and -2mV. Fi- 
nally we mention that the magneto-Peltier cooling via 
the bulk scattering can be also sensitive to the degree of 
energy relaxation, but discussion of the details is beyond 
the scope of the present paper. 

Since the thermopower-to-conductance ratios Si / Gi of 
the intermetallic interfaces studied up to now are smaller 
than the bulk values for thicker magnetic layers, for the 
material combinations considered above we do not ex- 
pect an increased cooling power by reducing the thick- 
ness of the nanopillars to the interface-dominated regime. 
The interface contributions are important for (classical) 
point contacts or pinholes in thick tunneling barriers, 
since Si can remain unmodified while Gi is strongly re- 
duced. Magnetic tunnel junctions are interesting sub- 
jects for magneto-Peltier studies since much higher S/G 
ratios can be expected. 

The spin Seebeck effect^ recently observed in a very 
long ferromagnetic metal appears to have a different 
origin than the conventional mechanisms of spin and 
heat diffusion. The observed thermoelectric spin sig- 
nal parametrized by a spin-Seebeck coefficient (S s = 
—2 nV / K at room temperature^) is much smaller than 
both the interface and bulk thermopowers considered 
above. We therefore do not expect that the spin See- 
beck effect would significantly modify our findings. 



VIII. SUMMARY AND CONCLUSIONS 

We studied the Peltier effect in nanoscale metallic mul- 
tilayer structures involving ferromagnets using a newly 
developed semiclassical theory of thermoelectric trans- 
port in magnetic heterostructures including spin relax- 
ations and the effects of electron interactions in limiting 
cases. The Peltier cooling/heating depends in general 
on the spin-degree of freedom as a function of spin and 
energy-dependent bulk and interface scattering. We pre- 
dict a magneto-Peltier effect in spin valves, i.e. a de- 
pendence of Peltier cooling on the relative alignment of 
the two magnetization directions, that can arise from the 
spin-polarization of thermopowers and is sensitive to the 
spin-flip scattering as well as strength of the inelastic 
collisions in the normal metal spacer. Similar behavior is 
found for the magneto-thermopower which might be even 
easier to observe in experiments (when thermoelectric 
voltage is measured rather than temperature). For fer- 



romagnetic layers with thickness of the order or smaller 
than the spin-flip diffusion length the magneto-Peltier 
effect should be observable in terms of magnetic-field- 
dependent resistance shifts in the R(I) characteristics 
i.e. the cooling-power. Estimates for the Peltier cooling 
based on our model and available parameters agree rela- 
tively well with experiments as well as numerical models 
in which the bulk scattering dominates. 
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APPENDIX A: PHONONS 

The Peltier effect in the presence of phonon heat con- 
duction and electron-phonon interactions can be mod- 
eled in linear response as follows: The net heat current 
flowing between the electron and phonon subsystems of 
the island for small temperature differences T™ — T™ -C 
Tp may be parametrized by the simple linear equation 
Q e - P = —Ke- P {Te — T™)r& For a phonon temperature 
drop of AT p across an interface, Q p — — k p AT p with k p 
the phonon thermal conductance of the junction. 34 The 
energy conservation laws then read: Q e i+Qe2+Qe- P = 
and Q p i + Q p2 — Q e - P = for the electron and phonon 
subsystems, respectively. The electron temperature in 
the node, Eq. ([S]), is then modified as follows 



K e i + K e2 + 7p(Kpl + K p2 ) 

and ATjf = jpAT™ where 7 P = K e - P / (k p i + K p2 + K e - P )- 
In the limit k p i(2) <C « e - P the Peltier cooling is reduced 
by the sum of the total thermal conductances Ku2) = 
K ei(2) + K P i(2)- The figure of merit SAT™ / AV is then 
further decreased by taking into account the contribution 
of the phonon heat conduction (Q p ^ 0) . 



APPENDIX B: THERMOELECTRIC SPIN 
DIFFUSION EQUATIONS 

In a diffusive magnetic metal in the steady state the 
Boltzmann transport equation in the relaxation time ap- 
proximation leads to the following spectral spin diffusion 
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equations for the local variation of the spin distribution 
functions /^^(e) for each spin a as 

V 2 / a (e)= /a(e) "{ 2 " a(C) . (Bl) 



where = y 1 D a r^ are the spin-dependent diffusion 

lengths. Under the detailed balance condition Af/ T J/ = 

A/"V r i/ the spectral spin diffusion equations can be 
rewritten as 

V 2 / s (6) = 4^, (B2) 

V 2 / c ( £ ) = -i 3 ^- (B3) 

where (l s f) 2 — (llj)~ 2 + (l'l,f)~ 2 and the charge and spin 
distribution functions, / c ( s ) = (/^ ± f* \ /2, have been 
introduced. In the strongly thermalized regime the spin 
diffusion equations can be expressed in terms of the spin 



chemical potentials /i-m) = (jtt c i Ms/2) and the electron 
temperature T. After inserting the linear expansions 

fs(e) « M„ (B4) 

V 2 / C (,)( £ ) - (-^) (vV W + ^V 2 t) , (B5) 

into the above diffusion equations we can integrate over 
energies by using the Sommerfcld approximation (see 
Eq. ([3])). We assume S 2 <C Co and disregard an en- 
ergy dependence of the spin diffusion length l s f (which 
is allowed when 2PS sf < (P' - P)S in which S sf = 
—eLoTd e lnl s f(e)\ ep ), but keep the energy dependence of 
the spin polarization P, recall Eq. (|2"2")) . One then arrives 
at the thermoelectric spin diffusion equations expressed 
in Eqs. (|39H41[) . Among the spin diffusion equations Eqs. 
([55)1 and (fJD]) are already well known^ Eq. (JUJ) repre- 
sents a spin-heat coupling for the electron spin diffusion 
in the presence of spin polarization of thermopower. 
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